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Abstract 

A microscopic theory for coarse graining diblock copolymers into dumbbells of interacting soft 
colloidal particles has been developed, based on the solution of liquid-state integral equations. The 
Ornstein-Zernike equation is solved to provide a mesoscopic description of the diblock copolymer 
system at the level of block centers of mass, and at the level of polymer centers of mass. Analytical 
forms of the total correlation functions for block-block, block-monomer, and center-of-mass pairs 
are obtained for a liquid of structurally symmetric diblock copolymers as a function of temperature, 
density, chain length, and chain composition. The theory correctly predicts thermodynamically- 
driven segregation of diblocks into microdomains as a function of temperature (c/ii parameter). 
The coarse-grained description contains contributions from density and concentration fluctuations, 
with the latter becoming dominant as temperature decreases. Numerical calculations for the block 
coarse-grained total correlation functions, as a function of the proximity of the system to its 
phase transition, are presented. Comparison with united atom molecular dynamics simulations 
are carried out in the athermal regime, where simulations and theory quantitatively agree with no 
need of adjustable parameters. 



PACS numbers: 61.20.Gy/61.25.Em/61.25.Hq/83.80.Sg 



I. INTRODUCTION 



One of the challenges in understanding the properties of polymeric materials stems from 
the necessity of developing theoretical approaches that can describe in a comprehensive 
manner properties observed at many different length scales. The presence of several length 
scales in which relevant phenomena takeplace leads to the complex nature of the liquid, 
rendering its treatment a difficult matter [if, 0,3]. Already in the description of the structure 
of homopolymer melts, two length scales need to be considered, which correspond to the 
monomer statistical segment length, cr, and the overall polymer dimension, i.e. its radius 
of gyration, Rg = cr(A^/6)^/^, where N is the total number of monomers in the chain. For 
diblock copolymers, the theoretical treatment is further complicated by the presence of a 
new length scale, which is intramolecular in character and corresponds to the size of a 
block. Diblock copolymers are macromolecules in which a homopolymer chain of = fN 
monomers of type A is chemically bound to a second homopolymer chain of different chemical 
structure containing Nb = {l — f)N monomers of type B, with = Na+Nb- The block size 
is defined by its respective radius of gyration such that, for example, the spatial dimension 
of the block composed of A-type monomers is given by RgA = a{NA/&Y^'^- 

Diblock copolymers are systems of great interest for their technological applications 
0, S]. Since the two blocks are chemically different, these experience a repulsive interaction 
that would encourage phase separation at low temperatures where entropy cannot balance 
enthalpic effects. However, the chemical bond existent between the two blocks prevents 
a complete separation of the two phases. As a consequence, at low temperatures block 
copolymer liquids undergo a microphase transition from disordered systems to ordered mi- 
crostructures of nanoscopic size, namely the microphase separation transition (MST). The 
length scale characterizing the size of the microphase is of the order of the block radius of 
gyration. 

With the purpose of developing the technology to produce micro-ordered structures of 
well-controlled size and shape, an understanding is required of the processes that drive the 
formation of micro-ordered phases under different thermodynamic conditions of temperature 
T and density p, as well as different chain composition /, monomer structure cr, and degree 
of polymerization A^. 

Computer simulations serve as an extremely powerful tool to investigate phenomena 
in complex fluids [g']. However, the limited power of present-day hardware does not allow 
for the simultaneous study of all length scales of interest. One way to overcome this prob- 
lem is through "multiscale modeling," where a set of simulations is performed at different 
levels of coarse-graining of the original system and in a subsequent step, information from 
different length scales is combined to provide the complete physical picture With such 
an approach, however, the challenge is not only to find the appropriate computational tech- 
nique for each length scale simulated, but also to know {i) the proper effective potential 
acting between coarse-grained units needed to carry out the simulations [s], [1, 10, ll[, and 
(m) the proper procedure for combining information from different scales of modeling once 
simulation data is acquired. 

First-principles theoretical models apt to coarse grain diblock copolymer liquids at 
different length scales of interest provide the potentials needed as an input in multiscale 
simulations, as well as the formal framework to combine information obtained from simula- 



tions of the liquid coarse grained at different length scales [12, [13!, [IJ, ll5|, ll6|. In a series 



of recent papers, we developed a coarse-graining approach that maps liquids of homopoly- 
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mer chains into liquids of soft colloidal particles, providing a formal analytical "transcale" 
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procedure [17|, [l8|, [l9(]. Each effective soft colloid is centered on the polymer center of mass 



and interacts through a Gaussian repulsive potential, which results in our formalism from 
the solution of Ornstein-Zernike integral equations. We later extended the same approach 



to describe the coarse-graining of homopolymer mixtures [17|, |20]. Computer simulations, 
performed by us, of coarse-grained polymer liquids and mixtures, where molecules interact 
by means of the derived effective pair potential, have been shown to reproduce quantitatively 
the structure and dynamics of the liquid at the center-of-mass level, while requiring con- 
siderably shorter computational time than that needed to perform united atom simulations 



In the present work, our approach is further developed to address the problem of mod- 
eling a melt of diblock copolymers as a liquid of interacting soft colloidal dumbbells. Each 
dumbbell represents one macromolecule composed of two effective soft colloidal particles, 
which in turn are sized according to the radius of gyration of each block and centered on 
center-of-mass coordinates of each block. Three different length scales are formally related, 
which correspond to coarse-graining the molecule at the monomer (the statistical segment 
length, a), block (the radius of gyration of block A, Rqa), and polymer (the polymer ra- 
dius of gyration, Rg) scales. In this way, our theory represents a minimal miramolecular 
mesoscopic model of polymeric liquid structures. 

There has bee n a g rowing interest in providing models for coarse-graining block copoly- 



mers chains. |21l . |22| . |23| | For example, building blocks of supermolecular structures, such as 
cellular membranes, could be modeled as self-assembling block copolymers chains. [2J] A re- 
cent paper proposes a model of coarse-graining for a symmetric diblock copolymer similar 
to ours, as the chain is modeled as two soft blobs, tethered by an entropic spring. j25| The 
blobs have equal size, and the coarse-grained total distribution functions are calculated nu- 
merically from a Monte Carlo simulation of diblock copolymers described at the monomer 
level. Monomers occupy the sites of a simple cubic lattice, with bond along the x-, y-, or 
z-directions. The two blocks individually are modeled as if they where in theta solvent, while 
the interaction between them is self avoiding. The numerical inversion procedure to derive 
the coarse-grained potential is performed in the athermal regime. As the authors point 
out in the paper, their model is "highly simplified", which proves the difficulty in treating 
intramolecular coarse-graining. The model, coupled with a reference interaction site model 
(RISM) and a random-phase approximation closure, predicts the mean-field clustering of 
diblock copolymers in a selective solvent. [26i] 

In this paper, we provide an analytical solution for the coarse-grained total distri- 
bution functions for a liquid of diblock copolymers represented as dumbbells of soft col- 
loidal particles. Our model differs from the one presented in Ref.([25|) in several ways. 
In our case, the size of the two "blobs" varies depending on the chain composition, /, 
degree of polymerization, N, and segment length, a. Moreover, repulsive interactions be- 
tween segments of different chemical nature are quantified by the interaction parameter, 
Xeff- Concentration-fluctuation stabilization enters through the polymer reference interac- 
tion site model (PRISM) theory for the monomer-level description. [271. 28l . 29 | and deviations 
from mean-field theory |30| are predicted by the coarse-grained approach as well. The two 
blocks follow Gaussian intramolecular statistics, which is a good approximation for copoly- 
mer melts, when each block has a degree of polymerization > 30, with a & A,B, and for 
the region in the phase diagram from the high-temperature to the weak segregation regime 
{XeffN << 10.5 for symmetric composition / = 0.5), where the system is isotropic. Numer- 
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ical mean field theory studies suggest coil stretching is not significant even below the order 
disorder transition until a strong segregation regime is entered, where XeffN > 100.|3ll. [si 



Analytical intermolecular total correlation functions between like and unlike coarse-grained 
blocks are predicted by our formalism as a function of chain composition, block size, density, 
temperature, as well as density- and concentration-fiuctuation screening lengths. 

One advantage of our approach is that it is analytical. Since a coarse-grained descrip- 
tion is obtained by performing statistical averages over local (small-scale) degrees of freedom, 
it translates the energy of the system into the free energy of the renormalized fiuid. jssi In this 
way, the obtained total correlation functions and related effective potentials, are functions of 
all characteristic physical parameters defining the system under consideration, such as tem- 
perature T, total site density p, and degree of polymerization N. For a diblock copolymer, 
the relevant parameters also include chain composition, /, and the interaction parameter 
Nxeff oc T~^. The parameter Xejj defines the proximity of the system to its order- disorder 
transition. For each set of parameters, total correlation functions and free energy change. 
As a result, a numerical solution of the coarse-grained description obtained from microscopic 
simulations requires performing a number of simulations equal to the number of combina- 
tions of those parameters, partially defeating the computational gains of a coarse-grained 
description. 

The material in this paper is organized in the following manner. We start in Section 
[IT] with a derivation from the Ornstein-Zernike equation of the total correlation functions of 
diblock copolymer liquids coarse-grained at the block level. Section IIIII provides a theoreti- 
cal description of diblock copolymer melts at the monomer level. In Section IIVI we present 
an analytically tractable solution in reciprocal space for a structurally symmetric diblock 
copolymer melt. The Fourier transform of the resulting expressions leads to analytical solu- 
tions of intermolecular total block-monomer and block-block pair correlation functions in real 
space, treated in Section |Vl Our analytical approach for coarse-graining the diblock copoly- 
mer liquid at the center-of-mass level and comparisons with the corresponding presentation 
of the homopolymer melt are discussed in Section IVI[ Finally, our theoretical development 
in the athermal limit is compared with united atom molecular dynamics simulations of a 
polyethylene melt in Section IVIH while temperature-dependent model calculations are pre- 
sented in Section IVIIII The paper concludes with a brief discussion and Appendices, where 
the auxiliary functions entering the exact solution of the total correlation functions, as well 
as the treatment of the block coarse-graining of a compositionally symmetric diblock, are 
presented. 



II. AN INTEGRAL EQUATION APPROACH TO COARSE-GRAIN DIBLOCK 
COPOLYMERS AT THE BLOCK LENGTH SCALE 

In this section, we derive the general expressions for the total pair correlation functions of the 
diblock copolymer liquid, coarse-grained at the level of the block length scale. The formalism 
is completely general and applies to any diblock copolymer system. The structure of a 
diblock copolymer liquid is characterized well by static correlation functions, which sample 
fiuctuations of monomer units in the fiuid. Given the position of monomer a belonging 
to block a G {A,B} comprising a chain j as r^, monomer fiuctuations at a specific wave 
vector k are represented by pi{k) = e^^''^°- . Since all chains are assumed to be equivalent, we 
henceforth discard the chain index j. Density fiuctuations for a generic monomer inside a 
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block a are defined as Pa(A;) = Z^^i e'^ '"'', while Pa{k)* defines its complex conjugate. The 
two-point density correlation functions, which describe the liquid structure, are given by the 
partial structure factor 
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which includes static correlations between monomers belonging to the same chain, Q'^^{k), 
i.e. the intramolecular static structure factor, and correlations between monomers belonging 
to different chains, H^^{k), i.e. the intermolecular structure factor. Since the liquid is 
spatially homogeneous and isotropic, the structure factors depend only on the modulus of 

— * 

the wave vector, |A;| = k. 

In our coarse-grained description, each block comprising the chain is mapped onto 
an effective particle, centered on the position of the block center-of-mass coordinate, Rba — 
Sa=i ^a- Fluctuations from the center of mass of block a are defined as Pbaik) = e**^'-^*"", 
and the partial structure factor becomes 



bb 



rbb 



(2) 



which requires knowledge of both intra- and intermolecular correlations. Block structure 
factors are normalized by the number of blocks in the chain, which is two for a diblock 
copolymer molecule. 

By analogy, block-monomer two-point correlation functions are given by 



S^{k) = -{pbaik) ■ ppikT) = + H^'Sik) , 



(3) 



which for compositionally asymmetric chains, f^j^ik) ^ f2p^{k), H^'p{k) ^ Hp^{k), and 
consequently -5^^ (A;) 7^ S^'^ik). 

Intra- and intermolecular structure factors are related through the Ornstein-Zernike 
equation, which has the general matrix formula 



u{k) ^ n{k)c{k)s{k) . 



(4) 



Here, C(A;) = fl^^{k) — S~^(A;) is the intermolecular direct pair correlation function matrix. 
In our description, the generalized Ornstein-Zernike equation includes contributions from 
"real" monomeric sites (m) and "auxiliary" sites positioned on the block center-of-mass site 
(6). 

At the block level description, matrices share similar arrangements. As an example we 
show the partial static structure factor, S(/c) = fl{k) + H(A;), which is defined as 



m 





g6m 




cmm 
'-'AA 
cmm 
^BA 


cmm 
'-^AB 
cmm 
'-'BB 


cbm 
'^AA 
cbm 
•^BA 


cbm 
'^AB 
cbm 
•^BB 


(^Qbrny^ 


gbb 




cbm 
•^AA 
cbm 
. i^AB 


cbm 
•^BA 
cbm 
'-'BB 


ebb 
'-'AA 
ebb 
^BA 


ebb 

•-'ab 

ebb 
'-'BB 



(5) 



where we omit the variable k to simplify the notation. In an analogous way we define the 
intramolecular structure factor matrix, which contains the correlation between real sites, 
^a/T — P'^Sia"' auxihary sites, f2^^^ = Ph^%, and the corresponding cross contributions. 
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n''^ = pbuj^^. Here, the block number density is pb = 2p/N, and cj^'U(^) = ^bb(^) = 1/2 
for a diblock copolymer. Furthermore, for blocks of different type, a ^ jS, we have that 
<^ik) = u^^^ik) and u%{k) = u;%{k), while u^^ik) + u'^ik). 

Consistent with its intramolecular counterpart, the matrix of total intermolecular pair 
correlation functions contains the correlation term for the real sites, HJ^^ = p^p^h'^g^, 
auxiliary sites H^^^ = pbaPb^h'^^p, and cross contributions H^*^ = pbapbph^. 
Here the number density of monomers A is pA = fp and = (1 — f)p, while the number 
density of blocks of type A (or B) is pbA = PbB = Pb/2 for a diblock copolymer. 

Finally, the intermolecular direct pair correlation function matrix includes the usual 
assumption that auxiliary sites are not directly correlated with either real or other auxil- 
iary sites. In this theoretical framework, we obtain the correlation in fluctuations of the 
intermolecular block-monomer function 



and the intermolecular block-block function 



1-1 



(6) 



(7) 



where Il'"^{k) = [U'^^k)]^ , since h''^{k) ^ hf'p^^k) whenever ^ with G {A,B}. 
The general block-block relation reads 



(8) 



where we used the property that both H'^'^(A;) and ri''"™'(fc) are symmetric matrices, along 
with the definition Qi^'^^k) = ^fl^'"^ {k)'j'^ . Eq.([7]) formally relates the total correlation func- 
tion for a coarse-grained diblock copolymer, represented as a dumbbell of two soft col- 
loidal particles, to the total correlation function and intramolecular structure factor of the 
monomer-level description. 

Eqs. ([6]) and ([8]) are further simplified using the fact that h^Xhi.^) ^^^baW 
(k) = Ii^baW- The block-monomer total correlation functions follows the general ex- 
pression 
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where 



and a(3 G AB^ while the block-block total correlation functions follows 
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Since no particular structure of the diblock copolymer has been assumed thus far, Eqs. 
(Q and (ITT]) are completely general and hold for any model of diblock copolymer chains, 
including diblock copolymers with asymmetric chain segments, as well as any general form 
of an interaction potential. From the knowledge of the pair correlation functions, obtained 
from the Fourier transform of Q(k) and H(A;), all thermodynamic properties of a polymer 
liquid can be formally derived [34]. 

In the following sections, we present an analytical solution for the intermolecular block- 
monomer and block-block correlation functions for a diblock copolymer liquid. We assume 
a structurally and interaction symmetric diblock with variable chain composition. The 
molecule is modeled as a Gaussian "thread" of infinite length and vanishing thickness. This 
model allows for an analytical solution of the total correlation functions in real and reciprocal 
spaces, as a function of the thermodynamic parameters of the system. 



III. MONOMER LEVEL DESCRIPTION OF DIBLOCK COPOLYMER LIQUIDS 

The coarse-graining formalism presented in Section [TTl is simplified when structurally and in- 
teraction symmetric diblock copolymers are investigated. For these model systems, segments 
of different chemical nature are assumed to have equivalent statistical lengths, o"a = cr^ = a, 
while the specific chemical nature of the block enters as an effective persistence length, 
through the block radius of gyration. Segments of like species interact through the poten- 
tials vaa ^ vbb, while unlike species repel each other through vab- At high temperatures, 
entropic effects dominate over enthalpic contributions, and block copolymer liquids resem- 
ble closely liquids of homopolymer molecules. As the temperature decreases, the effective 
repulsive potential Xeff = vaa + vbb — "^vab increases as Nxeff oc T~^, leading to the phase 
separation transition. This phase transition is characterized by a dramatic increase of the 
collective concentration fluctuation static structure factor, S'^{k*), at a specific length scale, 
k*. At the temperature of the phase transition, only certain fluctuations become anoma- 
lously large and the liquid segregates on a molecular length scale on the order of the overall 
size of the molecule, k* ~ Rg^- This remarkable property of copolymer liquids is due to the 
fact that, because of the connectivity between different blocks, even complete segregation 
cannot lead to macroscopic phase separation, as occurs in mixtures of two chemically differ- 
ent homopolymer melts. Since even at high temperatures, S'^{k) presents a peak due to the 
finite molecular size of the block copolymer chain, the peak position is largely independent 
of temperature. 

The first theoretical approach to describe the microphase separation transition was 



a mean-field theory developed by Leibler [30|. The theory is built on the expansion of 



the free-energy density of an ordered phase in powers of the order parameter, defined as the 
average deviation from the uniform distribution of monomers. The theory predicts a second- 
order phase transition for the symmetric mixture (/ = 0.5) and a first-order transition for 
asymmetric systems. Mean-field approaches are usually less accurate in the vicinity of the 
transition, where fluctuation corrections to the mean-field theory can change drastically 
the phase diagram. In the case of diblock copolymer melts, these corrections modify the 
predicted second-order phase transition for the symmetric diblock into a first-order phase 
transition. 

A fluctuation corrected mean-field approach was later derived by Fredrickson and 
Helfand [ssl] by implementing Brazovskii's Hartree approximation of Landau- Ginzburg field 
theory to treat diblock copolymer systems jSGj. The fluctuation corrected approach recov- 
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ers Leibler's results in the limit of infinite chain length. Fredrickson-Helfand predictions 
have been found to agree well with scattering experiments in the whole range of tempera- 
tures encompassing disordered to weakly ordered phases across the microphase separation. 
Since these approaches focus on the single-chain free energy and the hquid is assumed to be 
incompressible, fluctuation stabilization is mainly of entropic origin. 

Schweizer and coworkers developed an integral equation description of block copolymer 
melts that formally recovers the scaling behaviors obtained in field theories with only small 



differences [28|, 129|, |37|, 139|, |40]. In this case, however, the stabilization of the disordered 
state close to the MST is of enthalpic origin. Moreover, since the theory does not enforce 
incompressibility as a starting point in the treatment, the resulting structure factor contains 
contributions from both density and concentration fluctuations. This is consistent with 
the physical picture that in block copolymer melts far from their MST, the concentration 
fluctuation contribution is negligible while density fluctuations can still occur. 

In this work, we adopt an integral equation approach to describe the block copolymer 
structure at the monomer level, extending Schweizer's theory. This liquid-state description 
is largely compatible with the fluctuation-corrected mean-fleld approach, and has the ad- 
vantage of providing a theoretical framework that is formally consistent with the procedure 
presented in the previous section. In this way, the approach presented here builds on, and 
merges two well-developed theoretical flelds involving (z) the extension of integral equation 
approaches to treat complex molecular liquids [23, |4l[, and (ii) procedures to coarse grain 



macromolecular liquids [12|, Il5|, Il6|, ll7|, llSl, [19 



We focus on a structurally and interaction symmetric diblock copolymer. In the frame- 
work of an Ornstein-Zernike approach for this model system, repulsive mtermolecular in- 
teractions between A and B species, at the monomer level, are deflned by the direct pair 
correlation function containing hard-core intermolecular repulsive interactions between like 
species, Caa = Cbb = Cq, and a sum of repulsive hard-core and tail potentials for inter- 
molecular interactions between monomers of unlike species. Cab = Cba = Cq — Xejj/p- 
The effective interaction parameter, Xeff, controls the amplitude of microdomain scale con- 
centration fluctuations and increases as the system approaches its MST. Realistic diblock 
copolymer systems can be mapped onto this simplifled model, which has been extensively 
investigated in the past |4l|. In the theoretical coarse-graining approach presented here, the 
different chemical nature of the two blocks is accounted for by the difference in their radii 
of gyration, which is a function of the polymer local flexibility. This assumption is well jus- 
tifled in our approach since the monomeric structure is averaged out by the coarse-graining 
procedure, while local flexibility enters through the block radii of gyration. 

As a starting point in our derivation, we focus on the monomeric quantities which are 
input to our coarse-grained description for a diblock copolymer system, Eqs. and (fTTI) . 
As a flrst approximation, we assume that all monomers comprising a block are equivalent, 
so that each component in Eqs. and (ITT]) represents a site-averaged quantity. This 
is the conventional approximation adopted for treating analytically integral equations for 
polymeric liquids, and becomes correct when each block in the copolymer includes a number 
of monomers large enough to minimize chain end effects. The approximation is formally 
consistent with our coarse-graining theory, where physical quantities at the monomer level 
are averaged over the monomer distribution. 

The monomer total pair correlation function /^^^(r), with G {A,B}, in recip- 
rocal space is deflned as the difference between the total static structure factor and its 
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intramolecular contribution as 



p.p^Km = ^aT(^) - p</r(fc) , (13) 

where the total pair correlation function obeys the relationship 

We adopt here the thread model for the monomer-level description of the polymer 
chain. This model allows us to obtain analytical equations for our coarse-grained descrip- 
tion of the diblock copolymer system. In the thread model, the chain is treated as an infinite 
thread of vanishing thickness, with hard-core monomer diameter approaching zero, c? — 0, 
and a diverging segment number density in the chain, pmtra — ^ oo, while Pintrad? remains 
finite. The thread model yields a good description of properties at the length scale of 
including the presence of a correlation hole in the monomer pair correlation function. It 
cannot account for the local fine structure observed in the radial pair distribution function, 
g{r) = l + h{r), which is related to the presence of solvation shells due to monomer hard-core 
interactions. However, since our renormalized structures are characterized by a size compa- 
rable to the block domain, it gives a good representation of the coarse-grained structure for 
long chains where the block size is larger than the monomer diameter. 
Each of the site-averaged components of the total static structure factor contains contribu- 
tions from density S^{k) and concentration S'^{k) fluctuations, and can be expressed as a 
function of the density screening length as[39| 



Qmm(i,\ ^pi ^ ' 

^^^^^^ ^ + /s^{k) 

^^^^^^ ^ ^{i + k^e/fsm ■ ^^^^ 

The incompressible concentration structure factor is defined in Leibler's mean-field 
approach as 

= Fik) - 2Nxeff , (16) 

and diverges at the spinodal temperature, where Xejj = Xs = F{k*) / {2N) with the spinodal 
temperature defined as Nxs oc . The function F{k) = Nuj^^{k)/uj{k) depends on the 
intramolecular static structure factors ci;^(A;) with G {A,B} through the definition of 
uj{k) given by Eq. ffTOl) . and 

uZTik) = ooAAik) + 2u;^^{k) + a;^^(fc) • (17) 
To take into account the phase stabilization due to fluctuation effects, it is convenient to 



rewrite the incompressible structure factor in the following approximate form [43 



N/S%k) = F{k) - 2Nxeff = F{k) - 2Nxs + N/S^{k*) . (18) 

This expression takes into account the fact that when the spinodal condition is fulfilled, the 
inverse concentration contribution of the structure factor does not vanish: the disordered 
phase is still present and eventually the system undergoes a first-order phase transition. 
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With the purpose of obtaining an analytical expression for the coarse-grained sys- 
tem, we introduce the Otha-Kawasaki approximation |42j given by F{k) ^ (A/x + Bx — 
2[A5]V2) + 2Nxs, with X = eR^. Here, A = 3/[2p{l - ff] and B = l/[2/(l - /)]. 
At the peak position, defined as x = x* = [A/BY^"^ = {3/[/(l — f)]}^^"^, the contribution 
{A/x + Bx-2[ABY/^) ^ 0, thus recovering the spinodal condition of Xs = F{k*)/{2N). By 
inserting the Ohta-Kawasaki approximation into Eq. fllSp . the incompressible concentration 
fluctuation collective structure factor reduces to 

N/S'f'ik) ^ A/x + Bx- 2[ABY'^ + Ar/^^(F) , (19) 

which after expanding into partial fractions, leads to the following tractable analytical ex- 
pression for the concentration fluctuation contribution to the static structure factor jsoj. 



S^i.k) ~ 4h ^^-^ , (20) 



characterized by two length scales ^ = a[2BS'f>{k*) l?,]^'"^ and 6 = a{2BN /[U^ABf^ - 
3iV/5"^(A;*)]}^/2^ At the peak position, k = k* = {A/By/'^/Rg, the concentration structure 
factor behaves as S'^{k) = S'^{k*). In the small wave vector limit {k — >■ 0), it increases 
as S'^{k) oc k"^, in agreement with the observed scaling behavior for homopolymer liquids. 
Consistently, for the large wave vector limit {k ^ k*), it tends to zero as since S'^{k) = 
N/ (Bx), thereby recovering Leibler's scaling. The scaling with k at large and small wave 
vectors is also observed when homopolymer systems are investigated. This is a characteristic 
feature of block copolymer systems, signifying that at very large scales, as well as on very 
local scales, fluctuations are independent of the effective repulsion between monomers of 
unlike chemical nature. 

The solution of Eqs. iQ and ffTTl) relies on the definition of monomer-monomer and 
block-monomer intramolecular form factors. The monomer-monomer form factors are well 
described by the approximated Debye function. 



mm/7, N ^ 



k'^RgA 

2(1- f)NB 



k^Rg^ 



eRls-l + e-^'^'lB 



N 



^Asik) ^ ^ [k'RU - 1 + e-'=^"^-J [k'Rls - 1 + e-'^'^^- , (21) 



which can be conveniently simplified into their corresponding Pade approximants. [44 
In analogy with the center-of-mass monomer formalism [45], we approximate the block- 
monomer intramolecular structure factor in reciprocal space by the following Gaussian dis- 
tribution function 

a.^'^(fc) = (l/2)Ar^e-'='^-^/^ , (22) 

which includes the mean-squared radius of gyration, describing the squared average distance 
of a monomer of type (3 from the center of mass of an a-type block. 



R 



.2 



^ 1=1 



1 ^ 2 

-^(fft-i?J . (23) 
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In real space, uj^{r) is the generic intramolecular distribution function for any one of the 
Nf3 segments in a /3-type block with respect to the center of mass of an a-type block. Finally, 
we define the intramolecular total block-monomer structure factor as 

= ivTAik) + uj%{k) + u^Aik) + uj%{k) , (24) 

which in the A; ^ regime yields ^^7(0) = 1/2 [/ + (1 - /)] + 1/2 [/ + (1 - /)] N = N. 

While Eq. (1221) is a well-known expression when it applies to the distribution of 
monomers around the center-of-mass of an unperturbed homopolymer chain [i^, its ex- 
tension to diblock copolymers is novel. When tested against united atom simulation data 
(see Fig. [T]and the discussion of Section rVlT|) . these analytical expressions are fairly accurate 
for both symmetric and asymmetric diblock copolymers, while their simple Gaussian form 
allows us to derive analytical equations for the block total correlation functions in real space. 



IV. BLOCK COARSE-GRAINED DESCRIPTION IN RECIPROCAL SPACE, 
AND ISOTHERMAL COMPRESSIBILITY 



In the large-fc regime, which is of interest in block copolymer melts due to the finite size 
of the microphase separation, the ratios of intramolecular structure factors follow the rela- 
tions uj]^{k)/uj{k) ^ fcVV[12(l - /)], uj'g^{k)/uj{k) ^ A;VV[12/], and oj^]^(fc)/cj(fc) ^ 
l/[4iV/(l — /)]. By enforcing the approximation that uj^^{k) -C {i^aTI^)? '^rb which 
is justified by the fact that uj'^^{k)/uj{k) oc A^~^ and vanishes for long polymer chains, 
N ^ oo, the total block- monomer and block-block correlation functions simplify. Includ- 
ing these approximations into the monomer-monomer and block-monomer structure factors 
reduce Eqs. to the analytical general forms 



/^a"(fc)/2 



hp{k) + 
hp{k) - 



f-'u'^%k) - {1 - f)-'u;i^{k) 

uZT{k) 
f-'uj^ik)-il-f)-'uji^ik) 



(1 - f)Ah^{k) 
fAhHk) , 



(25) 



with a e {A,B}. 

Following the same procedure, the block-block total correlation functions in reciprocal 
space, Eqs. ( ITTjl . reduce to 



00 



tot 



[f-'ui^{k) - (1 - rr'ui^ik)] [r'ui^ik) - (i - fr'^i^ik)] 



UJ 



tot 



'{k)y 



f{l-f)Ah^ik) 
(26) 



where the density contribution hP{k) is identical to the monomer total correlation function 
for a homopolymer chain [27] = An^'p [^^/ (1 -|- ^pk'^) — ^'^/{l + ^c^^)] , and the concentra- 
tion fluctuation contribution at some thermal state point Nxeff oc T~^, having as a reference 
the athermal state Nxeff = 0, is derived from Eq. fl2U]) as Ah'f'{k) = hf^^^ ,{k) — h^ik) with 
h*{k) = [f{l-f)]-'p-'S^{k). 
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The total block-monomer intermolecular pair correlation function reads 



httik) = fhtA{k)/2 + (1 - f)h%{k)/2 + fh'S\{k)/2 + (1 - f)h%{k)/2 

^mm 

(27) 

where the contribution from concentration fluctuations rigorously vanishes, as is the case 
for the monomer level description of a diblock copolymer melt. The total block-block inter- 
molecular pair correlation function is given by 



+ 



f-' (uTAjk) + u'j.^^jk)) - (1 - {uj%ik) + g;^- (fc)) 



2 

Ah'f'ik) . (2J 



When compositionally asymmetric diblocks are investigated, the concentration fluctuation 
contribution to Eq. fl28|) does not vanish, but instead provides a small correction to the 
density fluctuation part. However, under athermal conditions or in the k —>■ limit, only 
density fluctuations are relevant since the concentration fluctuation contribution vanishes 
in a manner consistent with the monomer level description. It is worth noting that in the 
limit of a block approaching the size of the polymeric molecule, and in the limit of blocks 
of identical length (see Section [B] of the Appendix), Eq. fl251) recovers the homopolymer ex- 
pression for the molecular center-of-mass total pair correlation function, with concentration 
fluctuations strictly vanishing. 

As a test of our formalism, we present in Section IVIII a comparison of Eqs. (p5l) . 
(126|) and Eqs. (p7|) . (128|) against simulation data in the athermal regime. All equations 
show good agreement with simulations for both compositionally symmetric and asymmetric 
diblock copolymers (see also Figs. [2] and [3]), thus supporting the vahdity of our procedure. 

Finally, starting from Eqs. (1271) and (125]) . we calculate the isothermal compressibility of 
the system. Since the latter is a bulk property, it has to be independent of the level of coarse- 
graining adopted in the description of molecules in the liquid. The isothermal compressibility 
of the system coarse-grained at the block-monomer level, kt = S'l^ik 0)/{pkBT), is 
obtained from the matricial definition S'"^(A;) = f2'""(fc) + H''™(fc), after taking the — >■ 
limit and adimensionalizing the static structure factor. Each contribution is given by the 
relation Si^{k 0) = Ar^(e,VC')/2, which yields S'^^tik 0) = SZT{k ^ 0) = Nie,/^). 
In an analogous way, the compressibility of the system coarse-grained at the block-block 
level is calculated from the relation K-p = S^^{k — > 0)/{pkBT), and it is obtained from 
the matricial definition S^\k) = n''\k) + H''^(A;), where S^J'^ik ^ 0) = (^^/C^)/2. Since 
NS^t^^tik 0)/2 = S^t^{k ^ 0) = ^ 0) = N{il/il), this result is consistent with our 

prior findings obtained when coarse-graining homopolymer melts at the center-of-mass level, 
validating the feature that bulk properties, such as k-t, are independent of the fundamental 
unit (or frame of reference) chosen to represent the system. 

Reproducing the isothermal compressibility of the system, after performing a coarse- 
graining procedure, is an important test of the latter. Due to the nature of the coarse- 
graining procedure, effective coarse-grained potentials derived from pair distribution func- 
tions are softer than their microscopic counterparts. In fact, while real units, such as chain 
monomers, cannot physically superimpose, auxiliary sites can and the potential at contact is 
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finite. For this reason, the occurrence of small errors in the evaluation of the potential, which 
is often the case for numerically evaluated coarse-grained potentials, yields liquids that are 
too compressible. This shortcoming is eliminated in systems for which coarse-grained total 
correlation functions can be evaluated analytically, as it is in our case. 



V. ANALYTICAL BLOCK-LEVEL DESCRIPTION IN REAL SPACE 

A. Block- monomer total correlation function 

For a structurally and interaction symmetric diblock copolymer, the total pair correlation 
functions for block-monomer and block-block terms in real space can be expressed analyt- 
ically through a simple Fourier transform. The block-monomer expression separates into 
density and concentration fluctuation contributions as 

h':i{r) = hiy{r) + Ah^/{r), (29) 

with G {A,B}. The density fluctuation contribution is represented explicitly by the 
relations 

hT^'ir) = fFir, R,^,) + (1 - /)F(r, R,^,) = /.^'''(r) , (30) 

where for compactness, we introduce the auxihary function, Pir^R), deflned by Eq. (lAip 
of the Appendix. More specifically, Eq. (lAll) represents the density fluctuation contribution 
for one block, and is identical in form to the expression derived in our previous work for the 
center-of-mass-monomer total correlation function in homopolymer melts, coarse-grained at 
the polymer center-of-mass level Tsj . 



The concentration fluctuation contribution in real space is given by the relations 



fAh'Xfir) = Al^ir, R,aa) - AI%r, R,ab) = -(1 - /) A/.^"^"^( 

bm,()>f\ _ A p _ A TD ^ \ _ _ f \ h^m^'t' ( rr.\ 



r] 



(31) 



BA 



where each term is deflned as the difference in the response of the concentration fluctuation 
contribution between some thermal state {Nxejf oc T~^) and the reference athermal state 
(Nxeff = 0), as AJ-^ (r, R) = 4^^^. (r, R) - (r, R), with (r, R) deflned by Eq. of 
the Appendix. 

In the microscopic, small r/R regime, the concentration fluctuation contribution in 
Eq. (131]) reduces to 

A/^r, R) ^ ^fil - fr^e-'-'/^'^'^Nxeff • (32) 

This is the regime most relevant for block copolymer liquids approaching their phase tran- 
sition, since the microphase separation transition is characterized by segregation on spatial 
scales close in magnitude to the polymer radius of gyration. The temperature dependence 
enters Eq.( l3T]) through the Xejf parameter in Eq. (l32ll . evaluated at the reference athermal 
and thermal states. In this way, at high temperatures, i.e. far from the phase transition, 
density fluctuations are dominant over concentration fluctuations, and the total correlation 
function for diblock copolymer liquids recovers that of the homopolymer. 
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In proximity of the phase transition and for long polymeric chains, Eqs. (13T!) further 
simplify with the cross terms becoming negligible, while self terms yield the main contribu- 
tion, 

/A/.^-/(r) = -(1 - f)Ah'rAr) - a/^/^ (1 - fY^^e-'^'/^'^U.^Nx.ff, 

Vvr 2R,^^ ^^^^ 

V TT iKgBB 

Here, Rgap is the average distance of a monomer of type (3 from the center-of-mass of the 
block of type a, as defined in Eq. (12^ . A numerical study of these approximated expressions 
shows that neglecting terms with the most "cross" character is a reasonable approximation in 
real space, valid under different temperature limits and even when the system is asymmetric. 
A measure of the physical clustering with temperature is given by the parameter 

A/"^(r) = gtA{r) + g^^r) - g%{r) - g'^^ir) , (34) 

where g^{r) = 1 + h''^p{r). Since the number of monomers of type P included within a 
sphere of radius r' from the center-of-mass of block a is given by 



4vrp/. / r'g':^{r)dr , (35) 







with G {A,B}, Eq. flMl) represents a measure of the clustering of monomers around 
blocks of like species. The density fluctuation contribution to Eq. flMl) exactly cancels, while 
the concentration fluctuation contribution increases with decreasing temperature (increasing 
Nxeff)- At contact (r ^ 0), we obtain 

A^?'"(0) ^ yj/(l - f % {R-1a + R^bb) NXeff . (36) 

In the athermal limit block-monomer clustering due to concentration fluctuations vanishes, 
and Ag'^^^lO) = 0. At lower temperatures clustering of like species increases, with the 
leading factor being proportional to the ratio ^'p/Rga/3, which control the strength of finite- 
size coupling of microdomains {Rgap) and local (^^ correlations. 



B. Block-block total correlation function 



The block-block intermolecular total pair correlation functions can be solved in an analogous 
way to give in real space, 

h%{r) = h%%r) + Ah%^{r)^ (37) 

where the separation between density and concentration contributions is conserved. The 
density fluctuation contribution is given by 

= fJ' ir^Ro^ApA) + f{l ~ fV" {r^RaApB) 

+ /(I - f)JP (r, R^bpa) + (1 - ffJ' (r, Ro^b^b) , (38) 
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with G {A,B}, and J'^{r,R) defined by Eq. (lA2p of tlie Appendix. Tlie distance 
{Rap-ysf = [{Rga/B)'^ + {Rg'ys)]/'^^ wliere Rgap is tlie average distance of a monomer of type {3 
from the center-of-mass of the block of type a, as defined in Eq. (l23l) . We note that Eq. (lA2p 
was previously derived by us, in the context of coarse-graining a homopolymer liquid at the 
center-of-mass level IT], [l8|, [l9| . This term represents the total pair correlation function for 



a liquid of interacting soft colloidal particles, centered on the spatial position of the polymer 
center of mass. This simple analytical expression reproduces well data from united atom 
molecular dynamics simulations of polymer melts. 

The concentration fluctuation contribution is given by the general equation 

^h%\T) = AJ^ (r, R^apa) - ^«a/3b) 

- AJ^ (r, Ro^bpa) + AJ^ (r, RaBpB) , (39) 

where we define A J*^ (r, R) = jfixeg ~ "^o auxiliary function J^{r, R) by 

Eq. flASp of the Appendix. In the small r/Rg regime of interest here, the concentration fluc- 
tuation contribution simplifies, yielding for the generic contribution in Eq. fl5^ the following 
simplified expression 

A J^r, R) ^ f (1 - /)2y||e-3^V(4i?^) (^2 + 3^) Nxeff ■ (40) 

For long polymeric chains in general the cross statistical distances are larger than the 
self ones, Rga/s >> Rgaa, and Eqs. fl39|) simplify to 



V TT RgAA 



2 + 3 



(RgAA + RgBs)/'^ 




^h'^iir) - f (1 - fff-^/''^'^"''^^''''' I 2 + 3^ 1 Nxeff . (41) 

As with the block-monomer functions, numerical tests show that neglecting those "cross" 
terms is a reasonable approximation that holds under different temperature limits and even 
when the system is asymmetric. 

In analogy to the block- monomer development, to study the buildup of concentration 
fluctuations we define the parameter Ag'^^{r), which represents a measure of the physical 
clustering with temperature of blocks of like species, as 

^9'\r) = g%{r) + g%{r) - 2g%{r) . (42) 

The number of /?-type blocks included within a sphere of radius r' from the center-of-mass 
of block a, is given by 



(rO = 4np,p f r^g%{r)dr + , (43) 
Jo 
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with G {A,B}. Clustering due to concentration fluctuations increases with decreasing 
temperature, while density fluctuations provide a contribution constant with temperature, 
which is a consequence of the asymmetry in diblock composition and vanishes for compo- 
sitionally symmetric diblocks. The scaling with degree of polymerization of the function 
Ag^^''f'{r) oc ^/N / S{k*) depends on how far the system is from its microphase separation 
transition. At temperatures higher than the order- disorder temperature (T 3> Tqdt), we 
flnd that /^g^^^'t'{r) cx N'^/"^. At the transition temperature (T k. Tqdt), Ag''^''t'{r) oc N'^^^, 
while in the low temperature regime (T < Todt), A^*^''^(r) cx N'^^^. With respect to the 
monomer-block coarse-graining, A(7^''(r) contains a second term where the leading factor is 
proportional to the lengthscale ratio C,c/Rg- 



VI. COARSE-GRAINING OF DIBLOCK COPOLYMERS AT THE CENTER-OF- 
MASS LEVEL 

A. Reciprocal space representation 

In this section, we describe a diblock copolymer melt coarse-grained at the center-of-mass 
level. Information at this resolution is heavily averaged since the length scale of coarse- 
graining, Rg, is larger than the block size. However, it is still useful to consider this de- 
scription since it characterizes phenomena on the length scale of the polymer radius of 
gyration, and establishes a formal bridge of the theory presented here with previous ap- 
proaches to coarse grain homopolymer melts and their mixtures at the center-of-mass level 

To derive a coarse-graining procedure that maps block copolymer chains onto soft 
colloidal particles, centered on the coordinates of the polymer center of mass, the Ornstein- 
Zernike matricial relation is first specialized to include "auxiliary" center-of-mass sites. Here 
intramolecular structure factor matrix contains the matrix of real sites correlation defined 
before, as well as fl'^ = Pch^'^ with the number density of chain p^h = p/N, and ff^^ = pch- 
The matrix of the total pair correlation functions contains the elements H^™ = p'^hN^h'^ 
and H'^'^ = pl^h'^'^. The intermolecular direct correlation function matrix follows the usual 
assumption that there is neither a correlation between auxiliary sites nor with any other 
type of site, such that the only non- vanishing element is the contribution from C"^"^{k). 
Following analogous approximations and the analytical development of Section [HI we arrive 
to a representation of the relations cited above that rigorously decouples density and con- 
centration fiuctuations. This arrangement is given by the following set of expressions for the 
center-of-mass-monomer total correlation functions, 



^Z{k) 
WiTik) 



hP{k) + 
hP{k) - 



f-^uj--x{k) - {1 - fy'^'i.iky 

uZTik) 

7-Vl(fc)-(i-/)-V"^(fc) 

uZTik) 



(1 - f)/\h^{k) 

fAhHk) , 



(44) 



and for the center-of-mass total correlation function. 



h'^k) 



tot 



{k) 



ujZT{k) 



hP{k) + 



1, ,cm 



(k) 



f) 



-^uj'"^{k) 



uJ^Tik) 



/(I - f )Ah'^{k) , (45) 



16 



where we include the relation 



u^Zik) = u^^ik) + u^^ik) . (46) 
For the center-of-mass monomer intramolecular correlation function we start from the 



definition 45| uj^^^k) = Ne ^^^P^^ which leads to 



with 



R 



■gca 



— Y 



(47) 



(48) 



representing the radius of gyration involving A^^^ segments and the molecular center-of-mass 
coordinate, Rc- A justification for Eq. fHTj) can be found by performing the small- A; expansion 
of Eq. (j^6|) . In the athermal limit, Eqs. (jH]) and fHSl) correctly recover the homopoljTiier 



melt expressions [13, H 



In the case of a structurally and compositionally symmetric system, where io^iji) 



BB 



(k) and ^^(fc) 



B 
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lation for homopolymer melts [1 
[uj'''^{k)/u"'"'{k)]'^hP{k), as expected. 

The sum of domain-resolved contributions for the 



the equations further simplify recovering the known re- 
(k) ^ 



u;^™(fc)/u;™'"(A;)]/i^(A;) and h'^{k) 



intermolecular center-of-mass- 
monomer total correlation functions yields an expression analogous to the one obtained 
for polymer melts, where concentration fluctuations terms rigorously vanish, namely. 



hZik) 



UJ 



tot 



(k) 



h"''^{k) 



(49) 



Finally, the isothermal compressibility of the system at the present coarse-grained 
level, Kj- = S'^"^{k — ^ {)) / [pksT), can be obtained from the matricial definition S'^'"(/c) = 
f2'^™(A;) -|- H'^"^(/c), after taking the — > limit and adimensionalizing the static structure 
factor. The respective contributions are given by the relations 



S^^A^k 
S^"^{k 



0) = fme^/e) , 
0) = (1 - mic/e) 



(50) 



which yields after summing each contribution, that 5*^^ (A; ^ 0) = S'^^ik ^ 0) = A^(Cp/^c)' 
a result consistent with the expression for the liquid compressibility derived from coarse- 
graining homopolymer fluids at the center-of-mass level BHHIij. 



B. Real space representation 

An analogous representation is afforded in real space, where density (p) and concentration 
(0) fluctuation contributions separate as 
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with a G {A,B}. 

Since the integrands needed for the real space representation are identical in form to 
those presented in Section |Vl we simply cite the solution in terms of the respective functions. 
The density fluctuation contribution is represented explicitly by the relations 



h 



cm,p I 



fFir,R,^,) + il-f)Fir,R,^,: 



h 



cm,p I 
B \ 



(52) 



In the limit of ^ oo, the exact solution for the density fluctuation contribution, h'^"^'P{r) 



and h^'^'P{r)^ can be conveniently approximated |17l. Il8| by 



-3f2/2 



and 



39 
"16 



3 r 



^2 



9f 



"2 ~4 



-3f2/4 



(53) 



(54) 



where = ^p/Rg is the rescaled density fluctuation correlation length scale and r 
is the rescaled distance between intermolecular center-of-mass sites. 

The concentration fluctuation contribution denotes, as before, a difference between two 
thermodynamic conditions for the system under study, more specifically, between a thermal 
state {Nxeff oc T~^) and the athermal reference state {Nxejj = 0), yielding the expressions 



fAh 
Ah 



AI^{r,R,^,)-AI%r,R 



VcBy 



f)h 



B \ 



cc,(t)i 



AJ't'ir, Rg,A) - 2A J'^(r, Rg.AB) + A J^(r, R^.b) 



(55) 



where 27?^^^^ 



F?2 I d2 
^gcA "t" ^gcB- 



Numerical calculations of the mesoscopic correlations at the level of centers of mass 
show that h'^'^{r) is practically independent of temperature. This feature is consistent with 
the fact that intermolecular interactions between centers of mass of two block copolymers are 
unaffected by changes in concentration fluctuations, given that monomer correlations arising 
from two blocks are averaged out by the coarse-graining procedure. Moreover, the effect is 
intuitively explained by the fact that phase separation only occurs on the microscopic scale 
of Rg. As a consequence, Ah'^'^''^{r) k, 0. In contrast, the interaction between blocks, even 
for symmetric diblock copolymers, depends strongly on changes with temperature and the 
proximity of the system to the spinodal condition, as discussed in the previous sections. 



VII. NUMERICAL TEST OF THE COARSE-GRAINING PROCEDURE IN THE 
ATHERMAL LIMIT 

As a test of our coarse-graining expressions, we compare theoretical predictions with com- 
puter simulation data ^] of homopolymer melts in the athermal {Nxeff = 0) regime. 
Specifically, we use trajectories of united atom molecular dynamics (UA-MD) simulations 
of a polyethylene (PE) homopolymer melt composed of chains with degree of polymeriza- 

o _3 

tion N = 96, total site number density p = 0.0321 A , temperature T = 453 K, and 
Rg = 16.78 A [4^. Table [I lists the relevant length scales that enter into Eq. fl23|) . which 
are extracted from the simulation and used as an input to our calculations. By comparing 
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theory against simulations in the high temperature regime, where concentration fluctuations 
are not present, we can test the ability that our description has in capturing the effect of 
architectural asymmetry. We consider a diblock copolymer system where chain branches 
are of equal size (symmetric, / = 0.50), and where branches are of unequal size (asym- 
metric, / = 0.25), testing both intra- and intermolecular form factors at the block and 
center-of-mass levels. 

TABLE I: Length Scales for a Model Polyethylene Diblock Copolymer. 
Length [A] / = 0.50 / = 0.25 



RgAA 


10.86 


6.63 


RgAB 


27.80 


29.08 


RgBA 


27.79 


26.28 


RgBB 


10.88 


14.12 



While form factors in a diblock copolymer liquid at the monomer level have been 
extensively investigated, analytical expressions that represent well the structure factors at 
the block level are not known. As a starting point, we consider the block- monomer in- 
iramolecular form factors, which are assumed in this paper to follow Eq. (1221) . The latter 
is just a simple implementation of the well-known approximation for the center-of-mass 



monomer form factor in homopolymer melts [45|. In Fig. [H we test the Gaussian form 
of uj'^{k) per Eq. (122!) against simulation data for domain-resolved contributions. The top 
panel in Fig. [1] displays data for a compositionally symmetric diblock copolymer, for both self 
and cross contributions. The correlation between monomer and block center-of-mass sites 
decays faster in cross contributions, where the distance between the two species is larger 
than in the self contribution. The middle and bottom panels display the same comparison 
for data of a compositionally asymmetric diblock copolymer. We observe good agreement 
between the proposed expression, Eq. fl22l) . and simulation data for both symmetric and 
asymmetric diblocks. The Gaussian shape of the curve holds for any block-monomer form 
factor provided that the number of monomers in the block is sufficiently high, and the central 
limit theorem applies. The Gaussian form of the function allows for the analytical solution 
of the intermolecular block-block and block-monomer form factors. As a final check, the test 
of the total contribution, Eq. fl2^ . against simulations shows that the Gaussian form of the 
total intramolecular block-monomer structure factor also represents simulation data well. 

As a next step, we compare in Figs. [2] and H] our description of the mtermolecular 
structure factor at the block-monomer level with united atom molecular dynamics simulation 
data. We find that the agreement of analytical expressions, Eq. fl25l) and Eqs. (129|) to (1311) . 
with simulations is satisfactory in both real and reciprocal spaces. For the asymmetric 
system, / = 0.25, where the A-block is comprised of only 24 monomeric sites, the agreement 
between theory and simulation data tends to become rather qualitative whenever a site in 
the y4-block is involved, e.g. in h^XXi^) ^^"^ ^ab(^) largest discrepancy is localized near 
kRg^^ ~ 2 (Fig. [2]). It is interesting to note that while our approach predicts that h'^AA^k) and 
h^Xhi^) practically indistinguishable in the compositionally asymmetric diblock, data are 
sensitive to numerical errors due to finite-size effects. This discrepancy carries over to real 
space, where the data tends to be underestimated by the corresponding functions (Fig. Hj). 
However, the distinction between /i^(r) and h^Asir) is subtle, and our theoretical approach 
appears to provide a reasonable description both in real and reciprocal space. Moreover, 
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the h^shir) and h^BAi"^) terms are modeled rather well since the S-block is not affected by 
finite-size effects, and good agreement is also observed in reciprocal space. 

Next, we consider a comparison against simulation data of the theoretical block-block 
mtermolecular total correlation functions. In Fig. [31 we show the function in reciprocal 
space, hf'aisi^)' both compositionally symmetric (/ = 0.5) and asymmetric (/ = 0.25) 
diblocks. In both greement is found, within numerical error, between Eq. fl26p and 

simulation data. The top panel in the figure depicts the analytical solution involving both 
the Fade approximant of the intramolecular structure factors, as well as its Debye approx- 
imation, Eq. fl2T]) . The Debye approximation appears to give slightly better agreement 
with simulation than the Fade form. In turn, the Debye approximation yields only a minor 
improvement for the real space representations, as shown in Fig. [5l While the Debye approx- 
imation appears to model the data in slightly better fashion than the Fade approximant, the 
caveat in using it is that the reciprocal space representation must be numerically inverted, 
thereby defeating the purpose of obtaining an analytical solution. 

Shown in Fig.Elis the sum of contributions for the block-monomer total pair correlation 
functions. Upon inspection, it becomes evident that finite-size effects of the A-block in the 
asymmetric system cause a deviation from our theoretical predictions for r ~ 0.5 -R^^^, 
which is consistent with our prior findings for the domain-resolved functions. For larger 
length scales, however, the agreement is excellent. Discrepancies due to finite-size effects are 
not present in compositionally symmetric systems, which are modeled overall rather well by 
our theory. 

The sum represented by Eq. fl28p . which gives information of the liquid as a whole 
at the level of molecular blocks, is slightly sensitive to asymmetric features. However, 
the theoretical expressions are able to capture such small effects in both reciprocal and 
real spaces, as shown in Fig. [71 Specifically, blocks of different size induce a break in 
symmetry in the liquid, where packing is moderately favored at length scales smaller than the 
overall spatial extension of the molecule. While h^^{r) is identical for both symmetric and 
asymmetric cases, effects emerge at the block level that depend on differences in block size, 
a characteristic feature which enters in our development by the behavior of intramolecular 
uj'^{k) form factors. 

As a final test, we calculated the center-of-mass total pair correlation function given 
by Eq. f[32l) . The agreement between theory and simulations is good. Shown also is the 
discrepancy that arises when replacing uj^^{k) with Eq. together with the Gaussian 
approximation of the monomer-center-of-mass intramolecular structure factors, which results 
in a weak underestimate of h'^'^{r), as indicated in Fig. [3 

VIII. TEMPERATURE-DEPENDENT MODEL CALCULATIONS AND CLUS- 
TERING PHENOMENA 

In this Section we explore the temperature effects which give rise to concentration fluc- 
tuations in diblock copolymer systems. Our main goal here is to develop a qualitative 
understanding of the liquid behavior at the mesoscopic scale as the system evolves toward 
its microphase separation transition. To make contact with the calculations performed in 
the athermal regime, and presented in the previous section, we focus in our model calcu- 
lation on a "real" system, and we compute the cooling curves for the polyethylene system 
studied in Section [VTTl 

Input to our approach is the static structure factor, S^{k), described at the monomer- 
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level, as a function of the order parameter Nxeff- As discussed in Section llllt the static 
structure factor for a diblock copolymer liquid presents a peak that increases in intensity 
as the system approaches phase separation. At the temperature of the phase transition, 
only certain fluctuations become anomalously large, as the liquid segregates on the molec- 
ular scale of the molecular radius-of-gyration, k* ^ ^g^- Leibler's mean-field approach 
predicts a second order phase transition characterized by the divergence of the peak in the 
static structure factor, [S'^{k*)]~^ — 0. However, for experimental systems, where polymer 
chains are finite, the second order phase transition is suppressed by the onset of concentra- 
tion fluctuation stabilization, and a first order phase transition is observed over the entire 
composition range. 

At the monomer level, concentration fluctuation stabilization corrections are well ac- 
counted for by Brazovskii's correction to Landau approach, as described by Fredrickson and 
Helfand, as well as by the integral equation theory PRISM developed by Schweizer, Curro 
and coworkers. Both theories predict the same scaling with of the static structure factor, 
namely at high temperature (T >> Tqdt), random mixing applies and S'^{k*) oc N, while 
at the transition (T ^ Tqdt) they predict S^{k*)/N oc N^/^. The mean-field behavior is 
recovered in the limit of infinite chain length. The choice of the model used as an input, 
at the monomer level, is not crucial, however to preserve the consistency of our formalism, 
we correspondingly compute the cooling curve for the coarse-grained diblock in the frame- 
work of PRISM integral equation approach. The peak scattering intensity changes with 



temperature following the equation [371. |38 

t 



s'/\k*) [s'/\k*) -l] + [l-t] s{k*) -1 = 0, (56) 



where the form factor is normalized by its athermal value, s{k*) = S"^ (k*) / Sq (k*) , and the 
temperature, t = Thta,s/T, is rescaled with respect to the spinodal temperature obtained by 
applying the athermal, molecular PY closure. [28] Given the general relation Nx oc T~^, the 
temperature Thta,s is calculated in the framework of PRISM theory from the ratio Xs/ Xhta,s = 
(1 + Here a is the spatial range of the tail in the Yukawa potential governing AB 

interactions as fAB(^) = {a/r)ee~^^"', with e > being the interaction strength. For our 
calculations of polyethylene-like diblock, we set a = 0.5a, since this choice mimics the spatial 



range of the Lennard- Jones potential[37|, l38[, and = 0.346cr 



The other term entering Eq. that needs to be specified is 

iV 

where N = N{pa^Y is the Ginzburg parameter that controls the strength of the concentra- 
tion stabilization effect in S'^{k*). In addition, the parameter 

X* 

^ = c/(l-/)(2iVx.)V2' (^S) 

depends on = x*d'^F{x)/dx'^\x=x* which is of 0{1). Using the tabulated values from 
Ref. jssj, c = 1.102,1.278 and Nxs = 10.495,18.122 for / = 0.50,0.25, respectively. The 
calculation described so far is standard in PRISM theory and examples have been reported 
in several papers. 28, [2^, 37, H 
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The cooling curves, obtained from the self-consistent solution of Eg. (1561) for the system 
in this study, are presented in Fig. [H Note how concentration effects stabilize the response 
in S'^{k*)~^ for the finite-size PE system investigated, as evidenced by a leveling off upon 
further cooling. There is a subtle difference between the behavior for S'^{k*) for the / = 0.25 
and / = 0.50 cases, i.e. the two curves are indiscernible given the resolution in the figure. 
Also reported is the mean-field prediction, which shows divergence of the structure factor 
at the spinodal temperature. 

The system investigated is a diblock copolymer with fixed total number of monomers, 
= 96, identical segment lengths for the two blocks, cr^ = cb, and a repulsive Yukawa 
interaction between unlike monomers. The chain is partitioned, first as a compositionally 
symmetric diblock (/ = 0.5 and Na = Nb = 48) and then as a compositionally asymmetric 
diblock (/ = 0.25 and Na = "24: and Nb = 72). Input to our coarse-graining theory are 
the values of S^{k*) calculated for these two systems at Nxejj /Nxs £ {0.0,0.5,1.0,2.0}, 
as indicated in Fig. O These values sample our systems in a range of temperatures that 
include the athermal limit, Nxejj/Nxs = 0, the spinodal temperature, Nxeff /^Xs = 1, 
and the weak segregation limit down to (roughly) the ODT temperature, Nxejj /Nxs = 2, 



To study the effects of concentration fluctuations at the level of blocks, we focus on 
the physical clustering of like species as defined in Eq. (H2!) . In Fig. [TOl Ag''^{r) is shown 
for the symmetric and asymmetric cases. At athermal conditions concentration fluctuation 
contributions are absent. Repulsive interactions between unlike monomers are screened and 
entropic contributions to the free energy are dominant. The symmetric case exhibits no local 
clustering effects, and the system packs in an entirely random fashion. For the model diblock 
copolymer investigated in this study, where monomer bond lengths for the two blocks are 
equal {a a = ctb), the two blocks at high temperature are identical for a compositionally 
symmetric diblock, i.e. Ag^^{r) = for / = 0.5. 

For the asymmetric case, on the other hand, there is an emergence of entropic packing 
effects arising from the difference in block sizes, yielding a response in Ag^^{r) where AB 
contacts are favored ( A(7''^(0) < at high temperature for / ^ 0.5). 

This effect is quite apparent in the block frame of reference, which is sensitive to local 
microdomain scale correlations. Upon decreasing temperature, the formation of self con- 
tacts, AA and BB, becomes energetically favorable as the system approaches its microphase 
segregation transition (A(7''^(0) > 0). 

A shallow minimum develops with decreasing temperature, at a length scale corre- 
sponding to the size of the microdomain, r 1.5Rg for the symmetric case, because at the 
interface of the two domains the number of contacts between unlike species is higher than the 
number of self-contacts, i.e. Ag'^''{r) < 0. For compositionally asymmetric diblock copoly- 
mers, physical clustering occurs around the minority species, and the minimum is slightly 
shifted toward the small-r region. In both cases, the minimum is smooth and shallow, indi- 
cating that there is no a sharp transition at the interface between A and B domains, which 
is a characteristic feature of the weak segregation regime: fluctuations still partially disorder 
the liquid, while it becomes increasingly correlated approaching its phase transition. 

IX. CONCLUSION 

We have presented an analytical coarse-grained description for diblock copolymer melts, 
where the blocks in a polymer molecule are envisioned as two soft colloidal spheres connected 
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by an entropic spring. Corresponding domain-resolved intermolecular total pair correlation 
functions are formally derived from an integral equation approach by solving the Ornstein- 
Zernike equation, which is formalized as a matrix of monomer and block center-of-mass 
sites. The total pair correlation function describing the interactions occurring at the center 
of mass of the molecule is also presented. 

Analytical expressions for the total correlation functions of the coarse-grained diblock 
for a copolymer chain represented as a Gaussian thread of vanishing thickness, with an in- 
teraction symmetric potential acting between blocks of like and unlike chemical species, are 
derived. In the framework of this model, the analytical total correlation functions contain 
contributions from density and concentration fluctuations. The concentration fluctuation 
terms increase in intensity as the diblock melt approaches its microscopic separation tran- 
sition, however these do not diverge since finite-size fluctuations suppress the second-order 
phase transition. The contribution from concentration fluctuations drives the isotropic clus- 
tering of like species as the system approaches its phase transition. In the athcrmal regime, 
where density fluctuations are dominant, asymmetry in block size induces partial clustering 
of the domains. 

As a test of the coarse-graining approach, analytical expressions are compared with 
data obtained from a symmetric diblock copolymer melt in athermal conditions. Our the- 
oretical study shows that good agreement is attained in real and reciprocal spaces for the 
total intermolecular pair correlation functions at the block and center-of-mass levels. Com- 
parisons are made for a diblock copolymer melt composed of chains with equally-sized, as 
well as unequally-sized, chain branches. 

The present development corresponds to a signiflcant stride in presenting an analytical 
coarse-graining scheme for diblock copolymer melts. From our previous work, which has 
mainly focused on the mesoscopic treatment of polymers at the center-of-mass level, the 
results reported herein offer an intermediate level of coarse-graining that preserves some 
detailed intramolecular information to account for the physics proper of block copolymers. 
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APPENDIX A: AUXILIARY FUNCTIONS FOR REAL-SPACE REPRESEN- 
TATIONS 

1. Density Fluctuation Terms 

We collect here representations for the auxiliary functions used in the main text. The density 
fluctuation contribution arising from coupled frames of reference (i.e. those arising between 
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a mesoscopic level and the local monomer level) are represented by 
nr,i?) = -|(l-e.Ve3e«V(6«?) 



X 



e^/«^erfc|^ + ^Ue-^/^^erfc^ ^ 



(Al) 



Density fluctuations of self-character (i.e. contributions arising between the same mesoscopic 
level), are given by the relation 

R VS 



- e '^'■erfc 



+ 



(A2) 



2. Concentration Fluctuation Terms 



In the section, we summarize analj^ic expressions representing the contribution from con- 
centration fluctuations. Following the previous arrangement, fluctuations of mixed character 
are described by the expression 

/^r, R) = 3^5/(1 - /)| (I) ' e-^'-V^^^^) + /'^(r, k^, R) - I"f{r, k., R) (A3) 



with 



I'Hr, ^„ R) = T/(l - /)| (1 - e.kl) e«^4/a 



X 



e^^^erfc ( ^ + ^ 
a/6 V2R 



'■'=±erfc 









V2r) _ 



(A4) 



Analogously, the concentration fluctuation contribution of self character is represented by 
the relations 



J*(r,B)=/(l-/)|y'i|(|) .-»'V<4* 



2^^ 



r2 \ iV 
^~ZR2j ~-^^^~'^^5XF) ^ V /(I-/) 



J"^(r, A;+,i?) - J"^(r, A;_,i?) , 

(A5) 



with 



X 



e^^-erfc | ^ + e^^^erfc ^ "^^^ 



^/2> 2R 



2R 



(A6) 
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with 

k± = ^^'T^{^2r'. (A7) 

In these expressions the temperature dependence enters through N/S{k*) and the values 
of k± evaluated at the reference athermal and thermal states. The numerical evaluation of 
the complementary error function with complex arguments is more generally known in the 
context of Faddeeva's function in the field of optics, and poses no problem. However we 
note that, while the length scale associated with k± is complex, it is only a consequence of 
the factorization given by Eq. fl2U]) . In our expressions the imaginary components strictly 
vanish when taking into consideration the positive and negative branches of the functions. 



APPENDIX B: SIMPLIFIED FORMALISM TO COARSE GRAIN COMPOSI- 
TIONALLY SYMMETRIC DIBLOCK COPOLYMERS 



For a compositionally symmetric diblock copolymer (/ = 0.5), the general coarse-graining 
formalism presented in the main text becomes quite simple, since the equalities u^^lk) = 
uj'S'Sik), ^%{k) = w^l(fc), and uj''^{k) = uj%{k) apply. Also, we have that hJ'Xik) = 
h^'^{k), as well as h^AAi^) = ^bb(^)- By enforcing these rules, we obtain 



r^bm 



(k) 



h%{k) 



^tl{k) 
\k)\ 



ijj 



and 



Asik) 



tot 



^tTtjk) 
ik) 



-I 2 



L^tot 



^Zik) 



H 2 



h''{k) + 2 
hP{k)-2 

hf{k)+A 
hP{k)-4: 



u^Aik) - 


UJ%{k) 


ooTAik) - 


ik) 

^%{k) 




ik) 


'^TAik) - 


UJ%{k) 


, .mm 
^tot 

-u%{k) - 


ik) 

-u%{k) 



Ah'f'ik) 
Ah^{k) 



(Bl) 



-I 2 



Lo^rik) 



Ah'f'ik) , 
Ah'^ik). 



(B2) 



In real space, the block-monomer and block-block intermolecular total pair correlation func- 
tions separate into density and concentration fluctuation contributions. For the block- 
monomer contributions, density and concentration fluctuations become, respectively. 



^bm,p^ 



(B3) 



and 



Ah 



bm,q 
AA 



Ah' 



bm,4> I 



BB \' ) - 

-Ahf/{r) 



2Al\r,RgAA) 

4 

''BA 



2AI^{r,RgAB) 



-Ah'^lHr) 



(B4) 



Moreover, the density fluctuation contribution for the block-block correlation function is 
given by 



111 

- J^(r, Raaaa) + ^J^ir, Raaab) + ^ J^(r, Rabab) 



(B5) 
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while the corresponding concentration fluctuation contribution is 

Ahf^ir) = AJ'^ir, Raaaa) - 2A J'^(r, Raaab) + A J'^(r, Rabab) = -Ah'^i^r) . (B6) 

The functions IP{r, Rgap), J^{r, Rap^g), I'^{r, Rgap), and J'^{r, Ra^-fs) are deflned in Section 
|X]of the Appendix, with a, P,'j,6 G {A, B}. 
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FIGURE CAPTIONS: 



FIG. [T] Plot of uj'^{k). Shown are the Gaussian representations [lines] from theory 
and data from united atom molecular dynamics simulations correspondingly. The / = 0.50 
case is shown in Panel (a), whereas the / = 0.25 case is shown in Panels (b) and (c). Data 
is resolved into self- [circles] and cross-contributions [squares]. Panels (a) and (b), both 
display the self, uj'^{k), and the cross, u;^(A;), block- monomer contributions. In Panel (c), 
a;^(/c) is the self and uj'^{k) is the cross contribution. 

FIG. [2] Plot of h^f^{k). Shown are theoretical representations [lines] along with data 
from united atom molecular dynamics simulation [symbols]: AA [circles], AB [squares], BA 
[diamonds], and BB [triangles] contributions. Panel (a) is for / = 0.50, while Panel (b) is 
for / = 0.25. The dashed lines (which are indistinguishable in the plots) correspond to the 
solutions obtained from the Debye representation of co'^^(fc). 

FIG.ElPlot of h''^p{k). Shown are the theoretical representations [lines] along with data 
from united atom molecular dynamics simulations [symbols]: AA [circles], AB [squares], 
and BB [diamonds] contributions. Panel (a) is for / = 0.50, while Panel (b) is for / = 0.25. 
Panel (a) also shows the result from the Debye representation of uj^^{k) [dashed line]. 

FIG. S Plot of h'^{r). Data is arranged as in Fig. [3 

FIG. OPlot of Data is arranged as in Fig. [21 

FIG. E] Plot of [Panel (a)] and /i?™(r) [Panel (b)]. Lines are theoretical resuhs 

whereas symbols are data from united atom molecular dynamics simulations: shown are 
the / = 0.50 [circles] and / = 0.25 [squares] cases. The dashed line is as in Fig. [2l 

FIG. El Plot of hf^^{k) [Panel (a)] and hf^^ir) [Panel (b)]. Data is arranged as in Fig.l 

FIG. [H] Comparison of h'^'^{r) between theory and united atom simulation data for 
athermal conditions. Lines are theoretical results whereas the symbols are data from united 
atom molecular dynamics simulations. The representation obtained from the sum of uj'^^{k) 
terms is also shown [dot-dashed line]. 

FIG. M Plot of the cooling curves calculated from the integral equation approach, for 
a diblock copolymer system. Shown are results for / = 0.50 [solid hne], / = 0.25 [dashed 
line], and the mean-field behavior [dot-dashed line]. The points sampled for the model 
calculations are given by the circles. 

FIG. dU] Plot of ^g^^{r) as a function of the distance normalized by the 
polymer radius-of-gyration, for various temperatures. From bottom to top: 
Nxeff/Nxhta,s e {0.0,0.5,1.0,2.0}. Shown are the / = 0.50 [solid lines] and / = 0.25 
[dashed lines] cases. The arrows indicate the respective size of A-blocks. 
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